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Abstract 

We consider the rigid body quantization of Skyrmions with topological charges 1 to 8, as approx- 
imated by the rational map ansatz. Novel, general expressions for the elements of the inertia tensors, 
in terms of the approximating rational map, are presented and are used to determine the kinetic 
energy contribution to the total energy of the ground and excited states of the quantized Skyrmions. 
Our results are compared to the experimentally determined energy levels of the corresponding nuclei, 
and the energies and spins of a few as yet unobserved states are predicted. 



1 Introduction 

The Skyrme model [24] is a nonlinear effective theory of mesons, specifically pions. Its nonlinearity allows 
for the existence of topological soliton solutions, labelled by an integer-valued topological charge, B. A 
quantized Skyrmion of topological charge B is interpreted as a nucleus with baryon number B. 

The B = 1 Skyrmion was first quantized by Adkins, Nappi and Witten 0H]. They provided the first 
calibration of the Skyrme model, by fitting the model to the proton and delta masses. The toroidal B = 2 
Skyrmion was quantized in [HUE], and the energies corresponding to the ground state, representing the 
deuteron, and excited states were calculated. This analysis was extended in 20], to allow the two single 
Skyrmions to separate in the most attractive channel. This led to a more accurate determination of the 
mean charge radius, as the deuteron is rather loosely bound. 

The interpretation of the nuclei helium-3 and hydrogen-3 (triton) as quantized states of the B = 3 
Skyrmion was considered by Carson [8], and the spins and energies were calculated. This analysis was 
extended in [3] by a computation of the static electroweak properties of the quantized Skyrmion. 

In [28] , the B = 4 Skyrmion was semiclassically quantized, and the ground state (corresponding to 
the alpha particle) and first excited state were determined, and their energies calculated. The results, 
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though novel, involved consideration of selected vibrational modes as well as rigid body motion, and are 
difficult to generalize to higher baryon numbers. Here we will consider the B — 4 case from a different 
perspective, which may easily be generalized to higher baryon numbers, and enables us to compute the 
excitation energies of further excited states. 

Further results on the allowed spin and isospin states of quantized Skyrmions for B up to 8 and 
beyond have been obtained by Irwin [15], and taken further by Krusch [18) . However, in this work, there 
were no estimates of the energies of the states. 

It is not easy to assess the qualitative success of the Skyrme model just from the results for B < 4. 
The nuclei have the correct spin and isospin quantum numbers, but on the whole the ground states 
represent nuclei which are too small and too tightly bound. Nuclei with B — 2 or B = 3 have no excited 
states, experimentally, so Skyrmion excited states based on rigid body quantization are not meaningful, 
and one expects them to break up into individual nucleons if further degrees of freedom are included. 

There have been a number of developments which make it worthwhile to reassess these results on 
quantized Skyrmions, and it is also possible to extend them to the range of baryon numbers 1 < B < 8. 
First, it has been noted that a reparametrization of the Skyrme model is desirable to achieve a better 
fit to nuclear sizes and related quantities like moments of inertia [53]. The Skyrme length scale should 
be roughly doubled, and consequently the dimensionless pion mass parameter also doubled (to keep the 
physical pion mass fixed). Doubling the pion mass parameter has little effect on the qualitative character 
of classical Skyrmion solutions up to B = 7, but for B > 8 there is a clear difference [6]. The stable 
solutions are no longer the hollow polyhedra found earlier for B up to 22 and beyond, but instead more 
dense structures closer to what one expects for nuclei. In particular, for B a multiple of four, there are 
stable solutions which look like bound states of two or more of the cubically symmetric B — 4 Skyrmions 
[I]. We shall analyse below the quantum states of the B = 8 Skyrmion, which is made up of two B = 4 
cubes, and compare with the states of nuclei with B = 8, including beryllium-8. 

Another development is a better understanding of the quantization rules for Skyrmions, the so-called 
Finkclstein-Rubinstein (FR) constraints [TT], which encode the requirement that a quantized B = 1 
Skyrmion is a spin ^ fermion. The FR constraints combine the symmetry of a Skyrmion, for any value 
of B, with the topology of the Skyrme model, to constrain the spins and isospins of quantum states. 
Here the rational map ansatz comes in [T3] . This gives a separation of variables between the angular and 
radial dependence of the Skyrme field. True solutions do not exactly exhibit this separation, but they 
do so approximately. The ansatz gives a simple closed formula for the angular dependence of a Skyrme 
field, and rotational symmetries are easier to find than if one just has a numerical Skyrmion solution. 
The optimised rational map ansatz gives good approximations to true solutions up to B = 7 (and far 
beyond for zero or small pion mass). Even if it is a poor approximation, it can still be helpful in the 
numerical search for true solutions, and more importantly here, it is helpful in determining the effect of 
the FR constraints. Krusch has recently found a simple formula for determining the crucial signs that 
occur in the FR constraints [18 . This formula requires knowledge of the rational map approximating the 
Skyrmion. 

The rational map ansatz allows a further simplification, valid to the extent that a Skyrmion is well 
approximated by the ansatz. Kopeliovich noted that the moments of inertia (both rotational and isoro- 
tational) of a Skyrme field described by the ansatz are rather simpler than for a general Skyrme field 
[T7] . since the effect of a rotation is just to rotate the map, leaving the radial profile function invari- 
ant. We simplify Kopeliovich's formulae further, taking advantage of the complex analytic character of 
a rational map, and obtain formulae for the 36 components of the spin/isospin inertia tensor. These 
can be accurately evaluated, and it is easy to recognize if certain components vanish because of sym- 
metry. Using these moments of inertia, we estimate anew the energies of ground and excited states of 
quantized Skyrmions over the range of baryon numbers 1 < B < 4, and for the first time those in the 
range 5 < B < 8. The quantization is based on the established method of rigid-body quantization of 
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rotations and isospin rotations. Particularly interesting for us are the states of the B = 6 Skyrmion, 
since our reparametrization of the Skyrme model [23] was based on the mass and charge radius of the 
lithium-6 nucleus. Also interesting are the states for B = 8, since the double cube B = 8 Skyrmion has 
not previously been quantized. 

A problem for the Skyrme model that emerged in the work of Irwin |15j , is that the spin states of 
the B — 5 and B = 7 Skyrmions disagree with those of the corresponding nuclei in their ground states. 
It has been suggested more than once (see e.g. [2T]) that it might be appropriate, for these baryon 
numbers, to quantize a deformed Skyrmion with different symmetry. This would make sense, especially if 
the allowed spins were thereby reduced, making the spin energy smaller. The smaller spin energy might 
more than compensate the increased classical energy of the deformed Skyrmion. In this paper we are 
able to quantitatively assess this idea. For B = 7 it looks reasonable. A ground state with the correct 
spin | for the lithium-7/beryllium-7 isodoublet can be obtained, and the previously found spin | state 
interpreted as the observed, relatively low-lying second excited state. The spin ^ first excited state is 
still problematic, however. For B = 5 the situation is less satisfactory. 

In the next section we review the Skyrme model, and briefly describe the recent reparametrization 
of the Skyrme model using the lithium-6 nucleus [23) . Although this is very important, we show that a 
reparametrization alone cannot solve all the problems of the Skyrme model. In section 3 we describe the 
rational map ansatz for Skyrmions. Section 4 deals with the quantization of Skyrmions, which proceeds 
by parametrizing time-dependent solutions through collective coordinates. Here, we recall how the model 
is fermionically quantized by the imposition of FR constraints. In section 5 we present expressions for the 
inertia tensors which appear in the formula for the kinetic energy operator, in terms of the approximating 
rational map. Sections 6 to 13 deal with the energy levels of quantized Skyrmions of baryon numbers 1 
to 8 respectively. In section 14 we provide a conclusion. 



2 The Skyrme Model 

The Skyrme model is defined in terms of an SU(2)-valued scalar, the Skyrme field It is a low 

energy effective theory of QCD, becoming exact as the number of quark colours becomes large [2"0], 150]. 
We call the topological soliton solutions which it admits Skyrmions. 

The Lagrangian density is given by 

C = if Tr 9 ^ Ud>iU ~ 1 + ^ Tr [du.UU- 1 ,d v UU- 1 ] [d^UU-^d'UU- 1 ] + Tr (U - 1 2 ) , (1) 

where U(t,x) is the Skyrme field, F^ is the pion decay constant, e is a dimensionless parameter and m w 
is the pion mass. 

Using energy and length units of F^/4e and 2/eF^ respectively, we may express the Lagrangian as 
follows: 

L = y"|-iTr(i?^) + ^Tr([i? Al ,^][^,^])+m 2 Tr((7-l 2 )|d 3 x, (2) 

where we have introduced the su(2)-valued current i? M = (d fJl U)U~ 1 , and defined the dimensionless pion 
mass parameter m = 2m 7r /eF w . 

Field configurations of finite energy must satisfy the boundary condition U — > I2 as |x| — > 00. This 
compactifics R 3 to a 3-sphere of infinite size, and so topologically U : S 3 — * S 3 at a fixed time. Field 
configurations U therefore lie in topological sectors labelled by their topological degree 



B = J B Q (x)d 3 x, (3) 
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where 

B^x) = ^e^Trd'UU^d^U-'dPuU- 1 . (4) 

The degree B, which takes integer values, is identified with the baryon number. We refer to Bq as the 
baryon density. 

The kinetic part of the Lagrangian L is 

T = J {~Tr(R R )-~Tr([R i ,R ][R i) R })} d 3 x, (5) 

and this is quadratic in the time derivative of the Skyrme field. The rest of the Lagrangian ([2]) is (minus) 
the potential energy: 

E= J i-~Tr{R i R i )--^Tr{[R i ,R j }[R i ,R j })- m 2 Tr(E7 - 1 2 ) J d 3 x . (6) 

Static Skyrmion solutions can be obtained by solving the variational equations derived from E, or in 
practice by numerically minimising E in the sector with given B. 

The parameters e and F n can be fixed in a number of ways. It has been common practice to use the 
set of parameters given in [T] with the physical pion mass taken into account, specifically 

e = 4.84, F n = 108 MeV and m, = 138 MeV (which implies m = 0.528) . (7) 

In PP, the values of e and F n were tuned to reproduce the masses of the proton and the delta resonance. 
This parameter set was adjusted to optimise the predictions of the model in the B = 1 sector at the 
expense of the B = sector, which requires F v = 186 MeV. It is not, therefore, the optimal parameter 
set globally. 

In [23] , we proposed that in order for the Skyrme model to more closely model nuclear properties, a 
reparametrization would be desirable. We performed such a reparametrization by matching the model in 
the B = 6 sector with properties of the lithium-6 nucleus, obtaining 

e = 3.26, F n = 75.2 MeV and m„ = 138 MeV (which implies m = 1.125) . (8) 



Figure 1 shows graphs of nuclear masses and static Skyrmion masses per unit baryon number, using 
this new parameter set. The Skyrmion quantum energies are not included. We observe that the graphs 
intersect at B = 6, as expected. It is clear that it is not possible by a single parameter choice to correctly 
match nuclear and Skyrmion masses for all baryon numbers and this remains the case when the quantum 
spin and isospin energies are included. We believe that calibrating the model in the B = 6 sector is a 
promising way to describe the properties of nuclei with B > 4. For B = 1, 2, 3 the Skyrmion energies are 
now too high, and neither the nucleon mass nor delta resonance will be accurately fitted. 

For consistency, in the following sections we use the new parameter set ([8} throughout, making a few 
remarks about the old parameters in Appendix B. 
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Figure 1: Nuclear masses per unit baryon number (M/B) (solid), compared with static Skyrmion masses 
per unit baryon number (dotted). 



3 The Rational Map Ansatz 

We describe here the ansatz for Skyrme fields which uses rational maps between Riemann spheres to 
describe their angular behaviour [14] . This has been shown to give good approximations to several 
known Skyrmions, including all the minimal-energy solutions up to B = 7 (and much higher B when 
m = 0). The rational maps have exactly the same symmetries as the numerically known Skyrmions in 
almost all cases (B = 14 is an exception [5]). 

Via stereographic projection, the complex coordinate z encodes the conventional polar coordinates as 
z = tan(#/2)e J ^. Equivalently, the point zona sphere corresponds to the radial unit vector 



1 

1 + 1*1 

and inversely 



n z = — |_| 2 (z + z,i(f-z),l - |z| 2 ), (9) 

(n,)i+W a - (1Q) 
1 + (n,) 3 

The ansatz for the Skyrme field depends on a rational map R(z) = p(z)/q(z), where p and q are polyno- 
mials in z, and a radial profile function f(r). The target value R is associated with a point in the unit 
2-sphere of the Lie algebra of SU(2), given by the unit vector 

n R = T ^(R + R^(R-R),l-\R\ 2 ). (11) 

The ansatz is then 

U(r,z)=exp(if(r)n R{z yT), (12) 
where ri, T2 and T3 are Pauli matrices and f(r) satisfies /(0) = it and /(oo) = 0. 



Using this ansatz, the baryon number is given by 



B 



f/sin/Y /l 



2tt 2 V r J \l+\R\ 



dR 



dz 



2 



1% dz dz 

(i + M 2 ) 2 



r 2 dr , (13) 



and it can be shown that this is an integer equal to the degree of the rational map R. 
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The energy E for a field of form (|12D is 



E = 4tt 



2 f 2 + 2B sin 2 /(/ " 2 



■ 4 £ 

sin / 



in which X denotes the angular integral 



+ 2mV(l - cos/) dr 



(14) 



I = — 

47T 



1+1*1' 

l + \R\- 



dR 

dz 



2i dz dz 



(1 



(15) 



To minimise E one first minimises X over all maps of degree B. The profile function f(r) is then found 
by solving the second order ODE that is the Euler-Lagrange equation for the expression (fl"4"|) with B and 
X as fixed parameters. Given the profile function, the energy is determined by numerical integration. 
This gives the optimised rational map ansatz, and we denote the minimised energy by Mb- This is our 
estimate for the true Skyrmion mass, for baryon number B. To obtain a physical value to compare to a 
nuclear mass, one multiplies by the energy unit F^/Ae = 5.76 MeV, obtaining a classical Skyrmion mass 
in MeV. 



4 Quantizing the Skyrme Model 

Given a static Skyrmion Z7o( x ), there is generically a nine-parameter set of solutions, each with the same 
energy, obtained by acting with the Euclidean group and isorotations: 

U(x) = A 1 U (D(A 2 )(x - X))^ 1 (16) 

where Ai,A 2 are SU(2) matrices and A 2 is recast in the SO(3) form D(A 2 )ij — hTt(riA2TjA^ ). Semi- 
classical quantization is performed by promoting the collective coordinates A\, A 2 , X to dynamical 
degrees of freedom [7] . As we shall only be concerned with the computation of spin and isospin, we shall 
ignore the translational degrees of freedom X and quantize the solitons in their zero-momentum frame. 

Making the replacement U(x) — > U(x,t) = Ai(t)Uo(D(A 2 (t))x)Ai(t)~ 1 , and inserting this into the 
Skyrme Lagrangian, one obtains the kinetic contribution to the total energy 



T = ^ii.l ,,(/, - diWijbj + ^biVijbj , (17) 



where 

a,j = -il^YTjA^Ai , bj =iTiT ] A 2 A 2 1 . (18) 

b is the angular velocity in physical space, and a is the angular velocity in isospace. The inertia tensors 
Uij, Vij and Wij are given by: 

= -Jtt [T i T j + | [R^TiURk,^ d 3 x, (19) 

= - J eu m e jnp xix n Ti- ^R m R p + ^[R k , R m ][R k , Rp^j d 3 x, (20) 

Wij = J ejim xi Tr (r t R m + ~[R k ,Ti\[R k) rJ) d 3 x , (21) 

where R k — (dkUo)!!^ 1 is the right invariant su(2) current defined previously and Tj = | [rj, Uq] Uq 1 is 
also an su(2) current. The total energy, in terms of collective coordinates, is just T plus the constant 
M. b , the static mass of the Skyrmion. 
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We may write 

T= l -c T Wc, (22) 
where c T = (oi, a 2 , a 3 , b±, b 2 , 63), and the 6x6 symmetric matrix W is given by 

W =( v -^ T y J" (23) 

The momenta corresponding to bi and a; are the body-fixed spin and isospin angular momenta Li 
and K t 0: 

= U,>, • V.J'j- (24) 
A', = UijUj-Wijbj. (25) 

The usual space-fixed spin and isospin angular momenta Jj and Jj are related to the body-fixed momenta 

by 

Ji = -D{A 2 )j j L j , I i = -D{A 1 ) ij K j . (26) 
Defining H T = (Ki, K 2 , K3, L±, L 2 , £3), and using the relation H T = c T W, we find, provided det W^fl, 

T=^H T W~ 1 H. (27) 

We now promote the four sets of classical momenta introduced above to quantum operators, each 
individually satisfying the su(2) commutation relations. The Casimir invariants satisfy J 2 = L 2 and 
I 2 = K 2 . 

The basic FR constraints, which apply to any Skyrmion, are that physical quantum states |^) should 
satisfy 

^Trm-L^ = e 3wfa-K|^ = ) ( 2 g) 

for any unit vector n, which implies that for even B the spin and isospin are integral, and for odd B they 
are half-integral. There are further FR constraints on states if the Skyrmion has symmetries, and these 
are simple to determine if the Skyrmion is described by the rational map ansatz. A rational map, and 
hence the corresponding Skyrmion, has a rotational symmetry if it satisfies an equation of the form 

R(M 2 (z)) =M 1 (R(z)), (29) 

for some combination of SU(2) Mobius transformations M 2 and M\. M 2 corresponds to a rotation in 
physical space, and M\ to an isorotation. In general there will be a group S of such symmetries. We 
say that the map R is 5-symmetric if for each M 2 £ S, there exists an Mi such that ([29| holds. For 
consistency, pairs (M 2 , Mi) must have the same composition rule as in S, so R(M 2 M' 2 (z)) — MiM[(R(z)) . 
The map M 2 — > Mi is therefore a homomorphism. Note that it is not possible to construct such a map 
from Mi to M2. This is related to the fact that a Skyrmion may be invariant under a rotation alone, but 
cannot be invariant under an isorotation alone. 

Consider a rotation in physical space by an angle 9 2 about an axis 112, and an isorotation by an angle 
9i about an axis ni. We recall that under such a rotation, z transforms to M 2 (z), given by [TB] : 

= (cos I + i(n 2 ) 3 sin | ) z + ((n 2 ) 2 - »(n 2 )i) sin | 
(-(112)2 - «(n 2 )i) sin ^z + (cos ^ - i(n 2 ) 3 sin 4f) 

Similarly, under such an isorotation, R transforms to Mi(i?), given by: 

M (R) = (cos I + i(m) 3 sin %)R+ ((m) 2 - z(rn)i) sin % 
(-(ni) 2 -i(m)i)sin^ R+ (cos^--z(ni) 3 sinf ) ' 
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So given a specific symmetry (|29[) of a rational map, we use the above formulae to determine the corre- 
sponding angles and axes of rotation and isorotation. This is the data that is used in the conventional 
formulae describing the effect of rotations on quantized rigid bodies [19j . 

For #2 not an integer multiple of 27r, M 2 only leaves the points 

_ (n 2 )i+i(n 2 ) 2 _ -(n 2 )i-i(n 2 ) 2 

z n2 — — ana z_„ 2 — — — [6Z) 

l + (n 2 ) 3 l-(n 2 ) 3 

fixed. Similarly, Mi only leaves i?± ni fixed, where i?± ni are denned similarly. Therefore, for the symmetry 
1(25)1 to hold, we have 

#0-112 ) = Rai or #-n x ■ (33) 

Krusch showed that to correctly determine the FR constraint it is important to choose the direction of 
the axis ni so as to satisfy the base point conditior0 

#(z_ n2 ) = i2_ ni . (34) 

The symmetry (|25)l then leads to the following FR constraint on the wavefunction: 

e W a n,.L e ifl 1 n 1 .K|^ = XFfi |^) ; ( 35 ) 

a representation- independent statement, in which L and K are the body- fixed spin and isospin operators 
respectively, and the FR sign xfr = ±1. The combined rotation and isorotation corresponding to M 2 
and Mi respectively at most change the state by a sign factor, xfr- It was proved in [18] that the 
value of Xfr for a given symmetry of a rational map only depends on (9 2 and 0± , where the angles have 
unambiguous signs because of the base point condition (|34p , and is given by 

Xfr = (-1)^ , where N = ^-(B9 2 - ft) . (36) 

The FR signs Xfr form a 1-dimensional representation of the symmetry group S of the Skyrmion. 

A basis for the wavefunctions is given by | J, L 3 ) ® \I, K3), the tensor product of states of a rigid body 
in space and a rigid body in isospace. Here we suppress the additional labels J3 and ^3 which can take 
any values in the usual ranges allowed by J and I. J 3 is the physically meaningful projection of spin on 
the third space axis, and J 3 is the conventional third component of isospin. When we come to consider 
specific rational maps and the associated FR constraints, we seek low-energy states which are allowed by 
the FR constraints, and may represent the spin and isospin operators appearing on the left-hand side of 
(1531) by Wigner D-matrices, acting on the (2 J + 1) x (21 + l)-dimcnsional space of wavefunctions. 

Another advantage of the rational map ansatz is that it clearly illustrates any reflection symmetries 
of the Skyrmion, which enables one to determine the effect of the inversion operator V, which acts as 

P:U{x)^rf(-x), (37) 

and is baryon number preserving. For rational maps, the inversion x — * — x corresponds to z — > — 1/z, 
and the inversion U — > U* corresponds to R — * — 1/R. A rational map, and hence the corresponding 
Skyrmion, has a reflection symmetry if it satisfies an equation of the form 



-l/R{M 2 {z)) = M 1 {R{-l/z)). (38) 

In this case, V is equivalently given by the combination of rotation and isorotation corresponding to M 2 
and Mi occurring here. The parity of a quantum state is then the eigenvalue of the state when acted 



1 We note that this differs slightly from the condition given in I18| . This is because we are working with the inverse of 
the isospatial Mobius transformation which was considered there. 
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upon by the operator e * n 2' L e !9 i n r K derived from Mi and M\. There is, however, an ambiguity in 
the definition of V which was first explained in [15] . Given a candidate parity operator Vo for a given 
Skyrmion, we can also represent the operator by Vo times any element of the symmetry group of the 
classical solution. If the FR sign of a particular symmetry element is —1, then these two choices for V 
give different results. In particular, there is this problem for odd B: given a parity operator Vo, we can 
also represent the operator by Voe 2mnl ' , where n is any unit vector. As 2ir rotations have associated FR 
signs of —1 for odd B, we see that these two choices differ when acting on states. For each of the cases 
B = 1 to 8, we make particular choices for the parity operators, which we believe to be the most natural. 
We note that despite the ambiguity in the definition of V for a given Skyrmion, the relative parities of 
the Skyrmion's quantum states are fixed. 



5 Tensors of Inertia for Rational Map Skyrmions 

Kopeliovich [17] first presented general formulae for the inertia tensors of rational map Skyrmions. Writing 
Uq = exp(i/(r)n • r), these can be expressed as follows [T7] : 

Ua = 2 [ sin 2 / [{5ij - n inj )(l + f' 2 ) + sin 2 / 8^8^] d 3 x , (39) 
Vii=2 /sin 2 / 



(1 + f' 2 + sin 2 / d k n a d k n s ) (d m n r d m n r (r 2 ' b io - x l x ] ) - d l n r d j n r r 2 ) (40) 
- sin 2 / (d m n s d k n s d m n r d k n r (r 2 5 i:j - x l x ) - r 2 din r d k n r djn s d k n s ) d 3 x , 
Wij = 2 / eji m e lsp xin s sin 2 / [ (l + f' 2 ) d m n p + sin 2 / d k n r (d k n r d m n p - d m n r d k n p ) ] d 3 x . (41) 



These formulae for the inertia tensors assume that / depends only on r, and n depends only on the 
angular coordinates 9, <f>; further simplifications can be made if we assume that n depends just on a 
rational function R(z) as in (fTTj) . In order to obtain these simplified formulae, we find it helpful to write 
the R 3 metric and volume element in terms of r, z and z: 

ds 2 = dr 2 + r 2 d9 2 + r 2 sin 2 9 dtf = dr 2 + ^ = g a0 dx a dx f3 , (42) 

(1 + \z\ 2 ) 2 

,^ Ar 2 drdzdz , , 

dx = n . I 12^2 » 43 
(1 + \z\ 2 ) 2 

and to replace Cartesian derivatives with derivatives with respect to r, z and z. The products of com- 
mutators in (|19l2QI2ip may then be rewritten in these coordinates: 

[R k ,---}[R k ,---}^g rr [R r ,---}[R r ,---}+g zZ [R z ,---}[R z ,---}+g 2z [R- z ,---}[R z ,---], (44) 

where R z — {d z Uo)UQ l etc. We also have 

- itji m xiR m = (IjUo)!!^ 1 = HjR z - JljR z , (45) 

where 

- -H (i -* ,) s- (l -*>!)• (46) 



h - -^ 1 + ^ + ( 1 + ^lJ< (47) 
d _d 

h = Z—-Z-, (48) 
oz oz 
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so 



(49) 



We ultimately find that for the rational map ansatz, the tensors of inertia Uij, Vij and Wij can be 
expressed in the following form: 

*> - 2 / s " 2 ' H + ^ (lw I f |)> ' <»» 

where S = (£/, V, W) and the quantities . (which are given explicitly in Appendix A) are functions 
of the variables z and z only. In what follows, we use the above formula to numerically determine the 
elements of the inertia tensors, for a given rational map and profile function. The numerical values we 
obtain arc, of course, in Skyrme units. To convert to physical values we must multiply these by the 
mass scale and by the square of the length scale: (F„/4e) x (2/eF^) 2 — l/e 3 F 7T , obtaining quantities in 
inverse MeV. With the new parameter set ©, e 3 F„ = 2613 MeV. Although our optimal rational maps 
are familiar [13] , the profile functions have all been calculated anew using a shooting method^. In Fig. 2 
we plot the profile functions for B = 1 to 8, using the new dimensionless pion mass parameter m = 1.125. 




Figure 2: The profile functions f(r) for B = 1 to 8. B increases from left to right. 



6 B = l 

The rational map describing the single Skyrmion is given by R(z) — z, which is 0(3) symmetric. We 
find that the inertia tensors are each proportional to the unit matrix, satisfying Uij = Vij = Wij = XSij, 
where A is given by 

A = ^/rW/(l + /^^)dr. (51) 
Numerically we compute A = 45.1. 

2 due to Bernard M. A. G. Piette, University of Durham, UK 
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The FR constraints associated with spherical symmetry are 



where and n are arbitrary, leading to the following constraint on the space of physical states: 

(L + K)|tf)=0. (53) 

The "grand spin" M = L + K, or its components, appears quite frequently in what follows. The states 
satisfying (|5"3"f are linear combinations of the states |J, L 3 ) ® \I,K 3 ), whose grand spin is zero. These 
are of the form |J, I;M, M3) = \J, J;0, 0), in the standard notation for adding angular momenta, where 
J = L and I = K. So the spin J and isospin / have to have the same magnitude. In addition, J must 
be half-integral. 

The kinetic energy operator is given by 

T = —J 2 = —I 2 . (54) 
2A 2A V ' 

The eigenvalue of J 2 in states of spin J is J(J + 1), a standard result we will use frequently. Similarly I 2 
has eigenvalues /(/ + 1). For the lowest energy states, the nucleons with spin/isospin i, the spin energy 
in physical units is 

' ^e 3 ^ = 21.7MeV, (55) 



2(45.1)4 
and the total energy is 

E J=1 /2,i=i/2 =Mi + 21.7MeV = 986.2MeV + 21.7MeV = 1008 MeV. (56) 

This is not a bad fit to the nucleon mass 939 MeV. The spin energy is approximately one quarter of its 
value with the old parameters, but the higher classical Skyrmion mass makes the total energy too high. 
The spin/isospin | delta resonances come out with the too low energy 1095 MeV. Of course, using the 
old parameter set ([7]), the nucleon and delta masses are exactly right. 



As the rational map R(z) — z satisfies — l/R(z) — R(—l/z) 7 the parity operator V is naturally 
represented by the identity operator. We then find that each of the states described above has positive 
parity, in agreement with experiment. 



7 B = 2 



The symmetry of the B = 2 Skyrmion is -Doo/i , and the rational map which approximates this Skyrmion 
is R(z) = z 2 . The tensors of inertia Uij, Vij and W%j are all diagonal, with U\\ — U22, Vu = V22 
and Wn = W22 = 0. We also have that C/ 33 = 5W33 = IV33, relations which make the inertia tensor 
degenerate, a consequence of the axial symmetry. The degeneracy is resolved by imposing the following 
FR constraint on physical states: 

(£3 + 2^3)!*) =0. (57) 
The discrete symmetry R(l/z) = 1/R(z) leads to the FR constraint 

e** 1,1 ^** 1 !*) = -|*). (58) 

The ground state is then the J = 1, I = state |1,0) ® |0,0), which has the quantum numbers of 
the deuteron. The first excited state |0,0) ® 1 1 , 0} may be identified with the isovector 1 5o state of the 
two-nucleon system. 
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Using the expressions for the inertia tensors given in Appendix A, we find that numerically, 

U u = 96.58 , C/33 = 62.94 and V u = 160.61 . (59) 
The kinetic energy operator is given by [7J 

T = —J 2 + —I 2 -(— + - —) Kl . ( 60 ) 

2^ii 2U n \2U U V u W 33 J 3 V ; 

For the ground state, we find (with the conversion factor e^F^ implied from now on) 

Ej= x ,i=o =A1 2 + 16.3MeV= 1949.3 MeV + 16.3 MeV = 1966 MeV. (61) 

For the first excited state, we find 

£/=o,/=i = M2 + 27.1 MeV = 1976 MeV. (62) 

The experimentally determined mass of the deuteron is 1876 MeV, with the proton and neutron 
constituents only very weakly bound by 2 MeV. The state is marginally unbound, with a mass of 
1880 MeV As our energies (|6"Tj) and (f6"2"]) exceed the sum of the masses of a proton and a neutron, 
it would appear that we have predicted states that are unbound. However, when we compare Ej—± t i—Q 
and Ej = o. j = i to the sum of the masses of two quantized single Skyrmions with spin | (calculated in the 
previous section) , these states appear bound (with binding energies 50 MeV and 39 MeV respectively) . 

While the new parameters are clearly not ideal in the B = 2 sector, they predict results that are 
quantitatively quite accurate. The old parameters more strongly overestimate the binding energies of 
these two states [7J. Also, we calculate the excitation energy of the 1 Sq state to be 11 MeV relative to 
the deuteron, which is of the correct order of magnitude, and better than that obtained in [7] (35 MeV). 

To determine the parities of these two states we observe that the rational map R(z) = z 2 has the 
reflection symmetry —1/R(z) = — l/l), and so V = e mKs . Applying V to the allowed states, we find 
that both have positive parity, in agreement with experiment. 



8 5 = 3 



The tetrahedrally symmetric B — 3 Skyrmion was first quantized in [S]. Here we use the rational map 
ansatz to simplify the analysis. The Skyrmion is approximated using the map 



2 



z 3 — y/3iz 

The symmetry group is generated by two elements. These correspond to the following symmetries of the 
rational map: 

R{-z) = -R(z), (64) 
n / « + M iR(z) + l 

H^nJ = -iHM+i- (65) 

A 7r rotation about the X3-axis in space is equivalent to a tt isorotation about the 3-axis in isospace; 
and a 27r/3 rotation about the (x\ + xi + X3)-axis in space is equivalent to a 27r/3 isorotation about the 
(1 + 2 + 3)-axis in isospace. This leads to the FR constraints 

e i7TL3 e™ K3 \V) = I*), (66) 
e^( Ll+i2+L3 ) e ^( Kl+K2+ ^)|*) = I*). (67) 
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There is a spin |, isospin i (unnormalised) solution of these constraints, 

1 1\ 1 1\ 1 1\ 11 

17 2' 2/ 2' 2/ 2' 2/ 2' 2 



(68) 



This is the unique state with the same quantum numbers as the hydrogen-3/ helium-3 isodoublet of nuclei 

| and isospin |, 



in their ground states. The FR constraints also allow for two distinct states with spin I and isospin 3 



given by 



and 



I*) = 



I*) = 



3 3 
2' 2 



3 I 

2' _ 2 



3 1 
2' 2 



3 

2'' 



3 3 

2'~2 



3 3 

2'~2 



3 1 

2'~2 



3 
2' 



3 3 
2' 2 



3 1 

2'2 



3 

2' 



3 1 

2' 2 



3 1 

2' 2 



3 3 
2' 2 



3 3 
2' 2 



(69) 



(70) 



The first of these has the correct quantum numbers to allow for its interpretation as a nucleus in which 
one of the nucleons is excited to a delta isobar [TU] . 



The inertia tensors have been numerically determined for the rational map given above. They are all 
diagonal and proportional to the unit matrix: f/y = udij, Vij = vSij and Wi 
expected due to the irreducibility of the action of the tetrahedral group on R 3 . Numerically. 



wSij. This was to be 



u = 121.80 , v = 418.83 and w = -80.34 . 
The kinetic energy operator then takes the following form: 



1 



1 



T = - \(u - w)3 2 + (v- w)l 2 + wM 2 ] 

2 uv — w z J 



(71) 



(72) 



where M = L + K. Each of the three states (j68|69|70|) given above can be rewritten in terms of the 
basis states | J, I; . 

I 3 3. q o\ I 3 3. 



basis states |J, /; M, M 3 ): the first is proportional to | |, i; 0, 0), the second to 1 1, |; 0, 0) and the third to 



3, 2) — ||, |; 3, —2). They are thus eigenstates of M 2 with eigenvalues 0, and 12 respectively. The 
energies of the three states are then: 



^J=l/2, 7=1/2, M=o — 
Ej=3/2, i"=3/2,JW=0 = -A^3 ■ 



E 



J=3/2, 7=3/2, M=3 



3 u + v — 2w 
8 uv — w 2 
15 u + v — 2w 
8 uv — w 2 
3 5u + 5u + 6w 

8 uv — w 2 



-- M 3 + 15.4 MeV = 2895 MeV, (73) 
= /W 3 + 77.1MeV= 2957MeV, (74) 
M3 + 48.8 MeV = 2929 MeV . (75) 



These formulae are identical to those obtained in [S], although the numerical values of u, v and w 
are different because of the rational map approximation. The average mass of a helium-3 nucleus and a 
hydrogen-3 nucleus is 2809 MeV. Our ground state comes to within 4% of this value. However, our second 
state, with an excitation energy of 62 MeV, is rather too low in energy to have an NNA interpretation, 
which would require a mass splitting of approximately 300 MeV with the spin ^ ground state. Using the 
old parameter set, one obtains closer agreement with experiment. 

To determine the parities of these three states we observe that there is the reflection symmetry 
-1/R(iz) = iR(-l/z), and so V = e 1 ^ (^+^3). Applying V to the allowed states (|68I69I70|) . we find that 
they have parities +, + and — , respectively. We note that the helium-3 and hydrogen-3 ground states 
have positive parity. 
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9 B = 4 



The minimal-energy B — 4 Skyrmion has Oh symmetry and a cubic shape, and is described by the 
rational map 

= z 4 + 2V3iz 2 + 1 
z 4 - 2V3iz 2 + 1 ' 

This map has the generating symmetries 



= e^R(z), (78) 



R{iz) = 1/R(z), (77) 
+ 1 N 
-iz + 1. 

which lead to the FR constraints 

e i $S<. I 'i+ L '+ L ») e i 3 f-Ks\iJf) = (80) 

Seeking simultaneous solutions of these, we obtain the ground state |0,0) ® |0,0). There exists a spin 2, 
isospin 1 state given by 

(|2, 2) + V2i\2, 0) + \2, -2)) ® |1, 1) - (|2, 2) - y/2i\2, 0) + |2, -2)) ® |1, -1) , (81) 

and a spin 4, isospin state given by [T5] 

(V, 4) + ^|4, 0) + |4, -4) J ® |0, 0) . (82) 

The cubic symmetry excludes a spin 2, isospin state. 

The tensors of inertia are found to be diagonal, satisfying Un — U22, V%j — vSij and Wij = 0. 
Although the cubic group acts irreducibly on spatial M 3 , the associated isospin rotations are reducible, 
with the R 3 of isospace decomposing into a 2-dimensional and a 1-dimensional subspace. This is why the 
inertia tensor U has two independent diagonal entries, whereas V only has one, and why the cross term 
W vanishes. Numerically, 

U u = 142.84 , U 33 = 169.41 and v = 663.16 . (83) 
The kinetic energy operator is given by 

T = —J 2 + —I 2 + -(- )k 2 . (84) 

2v 2U n 2 \U 33 U n J 3 V ; 



For the spin 0, isospin ground state, the energy is simply the static mass of the Skyrmion, M4 = 
3679 MeV. Comparing this to the mass of the helium-4 nucleus, 3727 MeV, we see that our prediction 
comes to within 2% of the experimental value. The classical binding energy of the B = 4 Skyrmion is 
significantly larger than that of the B = 3 or B = 5 Skyrmion (see next section) . The mean charge radius 
of the quantized B — 4 Skyrmion was calculated using the new parameter set in to be 2.13 fm, which 
agrees reasonably well with the experimental value of 1.71 fm. Walhout [21] calculated this quantity using 
the old parameter set and taking into account a number of the vibrational modes, obtaining 1.58 fm. 

For the state (f5T|) with spin 2 and isospin 1, the energy is 

E J=2 , 7=1 = Ma + 28.7 MeV = 3679.0 MeV + 28.7 MeV = 3708 MeV . (85) 
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We note here that hydrogen-4, helium-4 and lithium-4 form an isospin triplet, whose lowest energy state 
has spin 2, and average excitation energy 23.7 MeV relative to the ground state of helium-4 [25], so here 
the Skyrmion picture works well. 



Finally, for the predicted spin 4, isospin state, we find 

Ej=a 7=0 = M A + 39.4 MeV = 3679.0 MeV + 39.4 MeV = 3718 MeV . 



(86) 



Such a state of helium-4 has not yet been experimentally observed. However, predictions for such a state 
with an excitation energy of 24.6 MeV have been made [T^l U3] • Our calculation suggests a slightly larger 
energy, in the range 30-40 MeV (allowing for the discrepancy between our calculation and the data for 
the isospin 1 state). The energy levels are summarized in Fig. 3. 



39.4MeV? 1=4*, 1=0 



28.7MeV_ J=2-J=l _28.7MeV_ J=2-J=l 

24.3MeV 1=2-. 1=1 



Hydrogen-4 



23.3MeV J=2-, 1=1 



28.7MeV 1=2", 1=1 



23.4MeV 1=2", 1= 



Lithium-4 



j=o + 1=0 

Helium-4 

Figure 3: Energy level diagram for the quantized B — 4 Skyrmion. Solid lines indicate experimentally 
observed states, while dashed lines indicate our predictions. 

To determine the parities of these three states we observe that the rational map (fTB"]) has the reflection 
symmetry — 1/R(z) = —R(—l/z), and so V = e l7rK3 . By acting with this operator on the physical states, 
we find that the spin 0, isospin state, and the spin 4, isospin state both have positive parity. On 
the other hand the spin 2, isospin 1 state has negative parity, and so we find no contradiction with 
experiment. 



10 B = 5 

Finding a quantized Skyrmion description of the ground and first excited states of the helium-5/lithium-5 
isodoublet, with spins | and |, has proved difficult. Physically, these states are not bound, and they 
may best be described as a cubic B = 4 Skyrmion loosely attracted to a single Skyrmion. 

It still remains to determine the symmetries and FR constraints that might give a lowest energy state 
of spin |. Here we explore in detail the idea floated in [5T], that one should consider variants of the 
rational map, and not just the one that optimises the classical Skyrmion energy. The minimal-energy 
B = 5 Skyrmion has Did symmetry, and it can be approximated by the rational map 

z(z 4 + ibz 2 + a) 
az A + ibz 2 + 1 



R( z ) = Z — 2 Z ^ , a = _3.07 , b = 3.94 . (87) 



The ground state obtained from this map [15] has spin | and isospin i, which is inconsistent with the 
observed spin | ground states of helium-5 and lithium-5. The Skyrmion has this shape up to a pion mass 
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in — 1 . However, if higher m is considered, the symmetry of the Skyrmion might change, and the above 
ground state become unstable. Therefore more symmetric solutions which apparently have higher energy 
might be the ones describing the true Skyrmion, and are worth investigating. 

When 6 = the rational map (|87p has D±h symmetry, and it acquires octahedral symmetry when 
in addition a = — 5. Octahedral symmetry has been previously considered [18] and it leads to a ground 
state with spin | and isospin i Let us therefore consider the D^-symmetric map (which could in fact 



be restricted to C4 symmetry) 
This map has the generating symmetries 

which lead to the FR constraints 



R \ z ) = J—-T> a^-b. 



R(iz) 
R(l/z) 



iR(z), 
1/R(z) , 



I*) 



e i7ril e i7rKl |*) = I*). 

Seeking simultaneous solutions, we obtain a ground state with J — | and isospin / = i, given by 



I*) 



3 




1 




3 




2' 


s) 


2 




2' 





1 1 

2'~2 



(89) 
(90) 



(91) 
(92) 



(93) 



This is the spin we are looking for. There are two excited states with J — § and I = \i most easily 
written in terms of | J, I; M, M3), where M = L + K: 



5 1-32 
2'2' 6 ' Z 



I i-3-2 

2'2'^' 



2'2' ' 



5I. 2 _2 

2'2' ' 



(94) 



with c± evaluated in Appendix C. The FR constraints allow for a further excited state with J = \ and 



-D4 symmetry implies that the tensors of inertia are diagonal, with U\\ — t/22, V\i — V22 and Wn 
W22, which leads to the expression for the kinetic energy operator 



1 



2 {U n Vn - W^] 



[U n (J 2 ~ L\) + V n (I 2 - Kj) + W n (M 2 J 2 I 2 2L 3 K 3 )] (95) 



(U 33 V 33 - Wi 3 ) 
The energy of the ground state is therefore 



[U 33 Lj + V 33 Kl + 2W 33 L 3 K 3 



Ej=3/2, 1=1/2 



3J7u + V n 



9U : 



33 T V33 



V 33 + &W : 



^UuVn-W?,) 8(U 33 V 33 - Wi 3 ) 



(96) 



The numerical value of the energy depends on the parameter a in the rational map, which has yet to be 
determined. 



We now argue that the D^h symmetry which we are considering is justified even if octahedral symmetry 
(a = —5) provides us with a slightly lower classical energy. The dependence of the classical energy on a 
is shown in Fig. 4, whereas the quantum energy is a strictly increasing function of a near a = —5 (sec 
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Fig. 5). Therefore the total energy achieves its minimum just below a — —5 (see Fig. 6), the quantum 
energy being much smaller than the classical one. Taking a = —5.0025 we find that 



U n = 203.41 , U 33 = 203.36 , V n = 1333.49 , V33 
Wn = -186.54, W33 = -186.61 . 



1332.96. 



(97) 



The static Skyrmion mass M5, for this value of a, is calculated to be 5101 MeV. The energies of the four 
states given above are then: 



^7=3/2, 7=1/2 

Ej=5/2, 7=1/2, c_ 
Ej=5/2, 7=1/2, c+ 



E 



,7=1/2, 7=3/2 



M 5 + 8.2 MeV = 5109 MeV , 
M 5 + 12.5MeV = 5114MeV, 
M 5 + 12.9 MeV = 5114 MeV, 
M 5 + 26.9 MeV = 5128 MeV . 



(98) 
(99) 
(100) 
(101) 



So, the achievement of the correct spin | for the ground state comes at a price. Firstly, the slightly 
excited J — \ state is not allowed by the FR constraints. Secondly, by comparing £7=3/2,7=1/2 to the 
average mass of the helium-5 and lithium-5 nuclei (4668 MeV), we see that our prediction is some 10% 
from the experimental value, whereas for the Z^d-symmetric Skyrmion there was an almost exact match 
to the helium-5/lithium-5 ground state energy. 



The Z?4/i-symmetric map (|88p satisfies —1/R(z) = R(—l/z). The parity operator could therefore be 
represented by the identity operator. However, we may also choose V = e 27rin ' L , where n is any unit 
vector. If we make this choice, then each of the states described above has negative parity. We note that 
the ground states of helium-5 and lithium-5 have negative parities. 
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Figure 4: Classical energy of the B = 5 Skyrmion as a function of a. 
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Figure 5: Quantum energy of the B = 5 Skyrmion as a function of a. 
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Figure 6: Total energy of the B = 5 Skyrmion as a function of a. 



11 B = 6 



The minimal-energy B — 6 Skyrmion has symmetry, and is well-approximated using the rational 
map 



R(z) = 



z + la 



z 2 (iaz 4 



1) 



This map has the generating symmetries 



which lead to the FR constraints 



R(iz) = 
R(l/z) = 



-R(z) , 
1/R(z) , 



(102) 



(103) 
(104) 



(105) 
(106) 



These constraints allow for the existence of states |1,0) <g> |0,0), |3,0) <g> |0,0), |0,0) <g> |1,0), |2,0) <g> 1 1 , 0) 
and |5,0) <gi |0,0). 

A numerical search over the parameter a in the rational map shows that the integral X is minimized 
at a = 0.16. However, it was suggested in [23] that allowing a slight deformation of the rational map 
would lead to more accurate predictions. In particular, it was found that by setting a = 0.1933, and 
using the new parameter set, one obtains a quantum quadrupole moment in agreement with experiment. 
In what follows, we set a = 0.1933. 

The inertia tensors have been computed for this rational map, and are each found to be diagonal, 
satisfying Uu = U22, V\\ = V21 and W\\ = W22 = 0. Numerically, 

[In = 215.84 , C/33 = 230.77 , V u = 1525.99 , V 33 = 1493.66 and W 33 = -105.45 . (107) 

The kinetic energy operator is given by: 

T = ^~ [ j2 - L l] + 2^ - X 3] + wu ,l_ w2 , [UssLj + V 33 Kl + 2W 33 L 3 K 3 ] . (108) 



2V n 



2(U 33 V 33 - ) 
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The static Skyrmion mass, Ma, is calculated to be 5601 MeV, which is precisely equal to the mass of the 
lithium-6 nucleus (the new parameter set was determined such that this would be the case - in [23j we 
estimated the spin energy for spin 1, isospin to be approximately lMeV, and then neglected this small 
quantity). The energy eigenvalues corresponding to the five states given above are then: 
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J=5, 1=0 
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Ma + 1-7 MeV = 5602 MeV , 
M 6 + 10.3 MeV = 5611 MeV , 



= M e 



12.1 MeV = 5613 MeV, 



— = Ma + 17.2 MeV = 5618 MeV , 



— =Ma + 25.7 MeV = 5626 MeV. 



(109) 
(110) 
(111) 
(112) 
(113) 



We may identify these with isospin states of lithium-6, and with states of the helium-6, lithium-6 and 
beryllium-6 nuclei, which together form an isospin triplet (see Fig. 7). The assumption in [23] that the 
spin kinetic energy of the state |1, 0) <8> |0, 0) is of order 1 MeV is clearly justified. 



15.5McV 1=2*,I=1 



10.4McV J=0*. 1=1 

Helium-6 



24.0MeV 1=5M=0 



15.5MeV J=2*,l=l 



10.4MeV J=Q*, 1=1 
8.6MeV 1=3M=<> 



IfoMeV J=2M=1 
10.4MeV 1=0*. 1=1 

Beryllium-6 



Lithium-6 



Figure 7: Energy level diagram for the quantized B = 6 Skyrmion. Energies are given relative to the 
spin 1, isospin ground state. 



Spin and isospin excitation energies relative to the lithium-6 ground state are experimentally known 
for these nuclei [Mj (see Fig. 8). The ground state of the lithium-6 nucleus is identified with the state 
|1,0) ® |0,0), and there is an excited state |3, 0) £g> |0, 0), with excitation energy 2.2 MeV. Lithium-6 has 
a further spin excited state with excitation energy 3.6 MeV; this is joined by the lowest energy states 
of the helium-6 and beryllium-6 nuclei, which relative to the ground state of lithium-6 have energies of 
4.1 MeV and 3.1 MeV, respectively, to form the spin isotriplet |0, 0) ® |1,0). Similarly |2,0) <g> |1,0) is 
identified with the spin 2 excited states of the isotriplet. A spin 5 excited state of lithium-6 has not been 
seen experimentally. We, however, predict the existence of such a state with excitation energy higher 
than those of the other states so far discussed. 

The splitting between the various spin and isospin states of the Skyrmion is clearly too large; the 
predicted quantum energies are roughly four times the experimental values. This may be connected to 
the fact that lithium-6 is an odd-odd nucleus. However, we have performed the same calculation using 
the old parameter set, and have found that this gives even wider gaps between the energy levels. The 
new parameter set is therefore not perfect, but is certainly an improvement. Furthermore, the ratios of 
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Figure 8: Energy level diagram for nuclei with B = 6. 



the relative excitation energies given by 

-Bj=o, i=i — Ej=i, i=o 
i?J=3, 7=0 — Ej=i, 1=0 

and 



= 1.2 (114) 



Ej =2 , I=1 -E J=1 , I=0 =h5 (n5) 
Sj=o, i=i — Ej=i t i =0 

correspond well to experimental data for these nuclei, for which the first ratio is 3.6/2.2 = 1.6 and the 
second is 5.4/3.6 = 1.5. 

To determine the parities of the states given above we firstly observe the reflection symmetry 



- l/R(e l -z) = -iR(-l/z). (116) 

The parity operator can therefore be represented as V = e % ~z 3 e~ % ? K3 . If we make this choice, then each 
of the states given above has positive parity, in agreement with experiment. 



12 B = 7 

Here, as for B = 5, quantizing the Skyrmion of lowest energy gives states with the wrong spins to match 
the nuclear data. The minimal-energy B — 7 Skyrmion has icosahedral symmetry, and is described by 
the rational map 

7z 5 + 1 
z 2 (z 5 - 7) 

This map leads to a ground state with J = |, I = i, a spin which appears experimentally as the second 
excited state of the lithium-7/beryllium-7 isospin doublet. Experimentally, the ground state has spin |. 



*(*) = * ■ 



There are many ways in which the icosahedral symmetry might be broken, allowing for the appearance 
ofaJ=|,/=i state in the spectrum. The most interesting possibility, in our opinion, is the breaking 
of the C3 symmetry, while preserving D5 symmetry. This leads to a ground state with J = | and I = |. 
So let us consider the Ds-symmetric map 

R ( Z )= 2f 5 +1 , ' ( 118 ) 

z z {z° — a) 
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where a = 7 restores the icosahedral symmetry. The generating symmetries of this map are 



R{e l ^z) = e-^R(z), 
R(-l/z) = -1/R(z), 



which lead to the FR constraints 
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and there exist two further excited states with J = |, 7 = 5, given by 

7 9 \ 1 1 \ 

I* 1 ) = 
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States with / = | are also allowed. In particular, there is one spin \ state 
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and two spin | states: 
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(119) 
(120) 
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The inertia tensors are found to be diagonal, with XJ\\ = U22, Vu = V22 and W\\— W22 — 0, leading 
to the kinetic energy operator: 



The energy of the ground state is given by 



[U 33 L 2 3 + V 33 K% + 2W 33 L 3 K 3 ] . (130) 



..31 9U 33 + V 33 - 6W 33 
Ej^/2, I= i/2 - M ^ + W^ + W^ + 8(U 33 V 33 - WK) ■ 



(131) 



The static Skyrmion mass, M7, is found to be close to 6328 MeV. The dependence of the classical and 
quantum energies on a are shown in Fig. 9 and Fig. 10 respectively. Looking for the value of a giving the 
minimum of the total energy, we obtain a = 7.002 (see Fig. 11). For this value of a, we find numerically 



U u = 246.27 , U 33 = 246.26 , Vn = 1873.03 , V 33 = 1872.76 , W 33 = 0.04 , 



(132) 
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Figure 9: Classical energy of the B = 7 Skyrmion as a function of a. 
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Figure 10: Quantum energy of the B = 7 Skyrmion as a function of a. 
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Figure 11: Total energy of the B = 7 Skyrmion as a function of a. 



£j=3/ 2 , 7=1/2 = M 7 + 6.6 MeV = 6335 MeV , (133) 
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to be compared to the average mass of the lithium-7 and beryllium-7 nuclei which is 6534 MeV. 
For the excited states we find that: 



#7=5/2,7=1/2 = M 7 + 10.1 MeV = 6338 MeV, (134) 

#5=7/2,7=1/2 = .M 7 + 15.0 MeV = 6343 MeV , (135) 

#}=7/2,i=i/2 = A^7 + 15.0 MeV = 6343 MeV , (136) 

#j=i/2, /=3/2 = M7 + 20.4 MeV = 6348 MeV , (137) 

#5=3/2,7=3/2 = M7 + 22.5 MeV = 6351 MeV , (138) 

#5=3/2.7=3/2 = ^7 + 22.5 MeV = 6351 MeV. (139) 

There are two main problems with the above spectrum. One is the absence of the J = 5 , / = 5 state, and 



the other is the appearance of the J = |, I = | state as the first excitation. We could try to overcome this 
problem by noticing that the first two excited states in the experimental lithium-7 and beryllium-7 spectra 
are in a sense anomalous: they have very low excitation energy, and the spin-energy correspondence is 
reversed. As was discussed in the introduction it is possible that such excitations cannot be described 
by our usual approach, and we need to allow for some vibrational modes or consider a Skyrmion of a 
different shape. Possibly, the states we find above could correspond to the ones lying above the lowest 
energy J = Z, I = i excited state. This interpretation fits rather well to the experimental data. The 
second problem is more difficult to tackle within this framework. The value of a being very close to 7 
leads to a configuration which is nearly C3-symmetric. This fact is reflected in the spectrum: we have two 
spin I and two isospin | states whose energies are almost indistinguishably close. This is not reflected 
in the experimental data. Let us therefore consider a smaller a and see if there is a better fit to the 
spectrum. Another advantage of this approach is that it helps to partially overcome the first problem as 
well. Indeed, by looking through a large range of a we find that at a — 2 the energies of the states given 
above are, in increasing order, 

X 7 + 6.3MeV, (140) 
X 7 + 9.3MeV, (141) 
A4 7 + 9.4MeV, (142) 
7W 7 + 13.7MeV, (143) 
7W 7 + 19.7MeV, (144) 
7W 7 + 19.9MeV, (145) 
jW 7 + 21.6 MeV. (146) 

The energy of the J = I, I = | state is now higher than the energy of one of the J = |, / — i states, 
and lower than that of the other, in agreement with experiment. This achievement, however, comes at a 
price as the classical energy is now some 10% higher than the experimental value. Figs. 12 and 13 are 
energy level diagrams for the quantized Ds-symmetric B = 7 Skyrmion, with a — 2, and for the B = 7 
nuclei, respectively. 

The symmetry breaking we have just considered is only one of the ways in which icosahedral symmetry 
might be broken. It is possible that a different breaking has to be considered in order to better understand 
the spectrum of the excited states, in particular the low- lying J = A, I = A state. It is also possible that 
with the increase of the pion mass the configuration will eventually break up into a B — 4 and a B = 3 
part, which is suggested by the very low energy for break up of lithium-7 into helium-4 plus a triton. 

The Ds-symmetric map (|118|) satisfies — 1/R(z) = R(— l/~z). As for B = 5, we find it favourable to 
choose V = e 27rm ' L , where n is any unit vector. If we make this choice, then each of the states described 
above has negative parity, in agreement with experiment. 
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Figure 12: Energy level diagram for the quantized B — 7 Skyrmion. A putative J = \ 
represented by dashed lines. 
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Figure 13: Energy level diagram for nuclei with B = 7. 



13 B = 8 



In this section we introduce some new ideas for estimating the moments of inertia of the B — 8 Skyrmion, 
and hence the excitation energies of the quantum states. It is believed that for our new parameter set, 
the minimal-energy classical solution resembles two touching B = 4 cubes (see Fig. 14) [4]. Here the 
rational map ansatz is not a good approximation, so our previous methods of calculation are no longer 
valid. Despite this, it is convenient to note that a field which is qualitatively of the right form, with the 
correct symmetries, can be obtained from a rational map, and this enables one to determine the allowed 
spin/isospin/parity states. There are also classical Skyrmion solutions which are well approximated by 
the rational map ansatz, and have only very slightly greater energy than the double cube. We consider 
these first. 
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For pion mass parameter between and approximately 1, the minimal-energy B — 8 Skyrmion has 
D§d symmetry, and is well-approximated by the rational map 



z 6 — ia 



which has the symmetries 



leading to the FR constraints 



^ = ^W^)' fl = ' 14 ' (147) 



R(e^z) = e^R{z), (148) 
R(l/z) = l/R{z), (149) 



e iwLl e ivKl \ty = |#). (151) 

The ground state is then determined to be |0, 0) (8> |0, 0) , and the first excited state is |2, 0) ® |0, 0) , in 
agreement with states of the beryllium-8 nucleus. However, this Skyrmion becomes unstable once the 
pion mass parameter exceeds 1. The true minimum is then described by two B = 4 cubes placed together, 
and as a first approximation to this in terms of a rational map we consider the O^-symmetric map (whose 
Wronskian vanishes on the 14 faces of a truncated octahedron): 

= z8 + A ^ z& - 10z4 + 4 ^ z2 + 1 fl 52 ) 
z 8 - 4V3z 6 - Wz 4 - As/Zz 2 + l' 



whose symmetries 



R(iz) = l/R(z), (153) 
1 -iz + l) l + V3R(z) 



lead to the FR constraints 



e^ L3 e ivKl \^) = |*), (155) 

Here the ground state is again |0, 0) ® |0, 0), but the |2, 0) <g> |0, 0) state is not allowed. The inertia tensors 
for this rational map are found to be diagonal, satisfying Uu — C/33, Vu = V22 = V33 and Wij = 0. 

However, the Oh symmetry is too strong for the description of two cubes, and has to be relaxed to 
D^h symmetry. Therefore we consider next 

. . z s + bz 6 - az 4 + bz 2 + 1 , lr . 

R ( z ) = -5 — r-« 1 — — 7 . ( 157 

w z 8 - bz 6 - az 4 - bz 2 + 1 y ' 

where a = 10 and b = 4\/3 restores the Oh symmetry. The rational map ansatz then gives a better 
approximation to the double cube Skyrmion, but only slightly because, for example, U = — 1 at the 
origin with the rational map ansatz, whereas for the true solution, U — — 1 at points near the cube 
centres. However, it has the right symmetry, and is good enough to determine the allowed spin/isospin 
states. The FR constraints are now 



*> = I*), (158) 
M*> = I*). (159) 
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which again allows a |2,0) ® |0, 0) state. The inertia tensors have the same symmetry properties as for 
the octahedral map, with the exceptions that XJ\\ ^ U33 and V33 7^ Vu = ^22- This leads to the kinetic 
energy operator: 

To determine the parities of states, we observe the rational map () 157(1 has the reflection symmetry 
— 1/R(z) = —l/R(—l/z). The parity operator can therefore be represented by V = e t7rI<2 . 

To progress, we now work directly with two cubic B = 4 Skyrmions separated along the Z3-axis, 
and find the moments of inertia of the resulting structure using the parallel axis theorem (ignoring the 
interaction of the cubes). The top cube is rotated by ? about the X3-axis relative to the standard 
orientation corresponding to (fT6")l . The bottom cube is rotated by — j about the X3-axis relative to the 
standard orientation. One difficulty here is in determining the separation of the cubes. The picture in 
Fig. 14 suggests that the separation is the value of r where the profile function becomes close to zero. 
From Fig. 2 we see that it is reasonable to take r — 1.8 leading to the separation in question being 
d = r I \f% =1.04 in dimensionless units. Then 



V±f= 8) = F 2 ( 2 S=8) = 2V/ 1 ( f- 4) + Md 2 = 2706 , (161) 

(B=8) _ 2V (b 
/ 33 — z ^33 



l/ 3 ( 3 B - 8) = 2V 3 f~ i] = 1326 . (162) 



where M. = 1277 (in dimensionless units) is the classical mass of two B = 4 Skyrmions. The isospin 
moments of inertia are simply given by 

r (B=8) _ TT (B=6) _ 07T (B=4) 



Ug=°> = ' =2Uii> =286, (163) 

t(B=S) _ tyjAB- 
y 33 — zt7 33 



U^~ 8) = 2C/^f" 4) = 339. (164) 



The equality of Un and U22, which we do not expect to be exactly satisfied by the true B = 8 solution, 
simplifies (|160p to 

T -wJ*-® + -mF-'® + & + &- (165) 

The ground state has quantum energy zero, so its total energy is simply the classical Skyrmion mass. 
The additional quantum energy of the spin 2, isospin state is 2.9 MeV, which is a very good match 
to the experimental value of 3 MeV 27J. There are a lot of further excited states, consistent with the 
FR constraints (|158p and (|159|) . whose wavefunctions, energies and parities are presented in Table 1. 
Figs. 15 and 16 are energy level diagrams for the quantized B — 8 Skyrmion, and for the B — 8 nuclei, 
respectively. We see a good agreement with experiment for positive parity states, and the appearance of 
some negative parity states which have not yet been observed experimentally. Of particular interest is 
the appearance of the J = 0, J = 1 negative parity state. If found, it could be a new ground state of 
the lithium-8 nucleus. The detection of the latter might be very difficult experimentally. We have also 
found quintets of / = 2 states. The lowest of these, with spin 0, have been detected experimentally with 
excitation energies very close to our predictions, and include the helium-8 and carbon-8 ground states. 

It remains worthwhile to find the optimal values of a and b in the rational map (|157p . Ideally, the 
classical mass should not be very far away from the experimental mass of the beryllium-8 ground state 
which is 7455 MeV, and the moments of inertia should be comparable with the ones we get from the 
double cube approach. The second condition is more difficult to achieve since the rational map is defined 
on a sphere, and cannot exactly reproduce a double cube configuration. Looking through a range of 
possible a and b values we find that the optimal map is given approximately by 

z 8 + 13^3 Z S _ 20z 4 + 13^ z 2 + 1 

R(z) = 2 -p 2 -p , 166 

Z 8 _ 1^3 Z 6 _ 2Qz4 _ 13^I Z 2 + !' 
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V n = V 22 


= 2901 , 


v 33 


= 2214, 


U u 


= 308, 


U22 


= 268, 


U 33 


= 283. 



leading to the following moments of inertia 

(167) 
(168) 
(169) 
(170) 
(171) 

We have recalculated the energies of the states in Table 1, using the kinetic energy operator (|160p and 
formulae for the energy levels of an asymmetrical top [19 . The energy of the first excited state is 

Ej= 2 , i=o=M$ + -^- = 7529 MeV + 2.7 MeV - 7531 MeV , (172) 

which is only slightly worse than the double cube approach. However, for further excited states the 
discrepancy in results increases, making the advantages of the double cube approach more evident. 




Figure 14: Baryon density isosurface for the numerically relaxed B = 8 Skyrmion with m ~ 1, resembling 
two touching B = 4 Skyrmions. 
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Figure 15: Energy level diagram for the quantized B = 8 Skyrmion, using the double cube approach. A 
putative J — 0~ isotriplet is represented by dashed lines. 
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Table 1: Energies and parities of the B 



allowed states. E^ c and E rrn are the quantum energies 



obtained using the double cube approach and rational map ansatz respectively. 
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Figure 16: Energy level diagram for nuclei with B = 8. 
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14 Conclusion 



The rational map ansatz simplifies the classification of the allowed spin and isospin states of quantized 
Skyrmions, and has enabled us to estimate their moments of inertia and energy spectra. The results 
are promising, and provide support for the interpretation of Skyrmions as nuclei. We have obtained the 
correct spin, parity and isospin quantum numbers for the ground states and various excited states in most 
cases, and the quantum energies of excited states are reasonably close to the experimental values. We 
have also been able to predict some excited states that have not yet been observed. The new parameter 
set for the Skyrme model, with which we have been working throughout, has provided better results than 
the traditional parameter set for the larger values of B. We have also put into effect a new approach 
for some Skyrmions of odd baryon number, in particular for B = 7. By deforming the highly symmetric 
minimal-energy Skyrmion, we have been able to reproduce the spins of the experimental ground state 
and several excited states. We have given the first estimates of the energies of quantum states based on 
the double cube B = 8 Skyrmion, and similar methods should be applicable to the multi-cube solutions 
for B = 12, 16 and beyond, presented in [I]. 

The calculations presented here are subject to a number of limitations. Firstly, we consider the 
semiclassical quantization, in which only the collective coordinates for rotations and isospin rotations 
are considered. A more accurate procedure would have to take into account further degrees of freedom, 
which we refer to as vibrational modes. Our current understanding is that the Skyrme model provides 
a description of nuclear physics in which nuclcons are partially merged, and their orientations in space 
and isospace are highly correlated. In a sense, this is the opposite of a naive shell model, in which 
nucleons move in a potential, and are to first approximation uncorrelated. A more realistic model would 
possibly lie somewhere between these two extremes. Allowing the individual Skyrmions, or subclusters 
of Skyrmions, to move relative to each other, and performing a quantization of these degrees of freedom, 
would be a significant refinement to our approach. In so doing, some missing low-lying experimentally 
observed states of nuclei may appear. These include the low- lying excited states with J = \ and I = \ 
that are present for B = 5 and B — 7. 

The work here should be taken further by working with the exact Skyrmion solutions, and not just 
the rational map approximation to these solutions. Classical energies and moments of inertia will change, 
though we hope not drastically. Further investigation of the effect of varying the dimensionless pion mass 
parameter is also warranted. The length scale of the Skyrmions is quite sensitive to this. Possibly, an 
increased parameter will create an instability in the Skyrmions with B = 5 or B = 7, thereby justifying 
our arguments for changing the symmetries. 
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A Appendix: Inertia Tensors 



The tensors of inertia for rational map Skyrmions may be expressed in the form: 



Ey = 2 J sin 2 /- 



C'y 



\R\ 2 ) 2 



2 sin 2 / ( l + \z\> 
J r 2 U + |i?| 2 



dR 



dz 



where X = (U, V, W) and the quantities Cu tj are given by 







= \l-R 2 \ 2 
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Cu 22 


= \l + R 2 \ 2 
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= 4|i?| 2 , 




c Ul2 


= c U21 


= -2 3?i? 2 , 




c Ul3 


= c U31 


= 2(|i?| 2 - 


1) MR, 


c u 23 


= Cu 32 


= 2(|i?| 2 - 





the quantities . are given by 



C Vll = 



2 1 2 



dR 



Cv 22 
Cv 33 
C Vl2 = C V21 
C Vl3 = C Vsi 
Cv 23 = Cv 32 



\l + z 2 



dz 
dR 



4\z 



dR 



dz 
2 3z 2 



dz 
2 



dR 



dz 



2M(\z\ 2 z-z) 
23?(|z| 2 z + z) 



dR 



dz 
dR 



dz 



d 6 x. 



(173) 



(174) 
(175) 
(176) 
(177) 
(178) 
(179) 



(180) 
(181) 
(182) 
(183) 
(184) 
(185) 
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and finally, the quantities Cw i6 are given by 



O,,, »M !J -z 2 )(l-R 2 )^\ , (186) 

<V,. >M (l + * a )(l+5 a )^) , (187) 

r,i ... -i )i x ( £*^J , (188) 

r,,,, = -3 ( (i + z 2 )(1 _ j r2 ) ^ ; (189) 

r U| , = -2 l^(l-fl 2 )^J , (190) 

r,,,, = . i :il • /. ,J i^) , (191) 

fu,, = 3 ( ( I - 2 2 )(1 + j r 2) ^ ^ (192) 

O, , = -2 1?(i2(l-* 8 )^ , (193) 

O, , - 2o j //(l + z 2 )^ . (194) 



B Appendix: Old Parameters 

Here we collect some data on moments of inertia, in Skyrme units, calculated with the dimensionless 
pion mass parameter m = 0.528 that emerges from the calibration of [T]. The following results are novel, 
as they were obtained using the rational map ansatz and the formulae in Appendix A, and extend from 
B = 1 up to B = 4. For B = 1 the rational map ansatz is exact, so our result should agree with that of 
PQ, and indeed it does. For B — 2,3 our results can be compared with the moments of inertia calculated 
from the exact Skyrmion solutions (with the same m) as given by [71 [5] • This allows us to investigate the 
accuracy of the rational map ansatz for these Skyrmions. 

The notation is as in the Sections 6 to 9 above. For B = 1 

A = 62.85. (195) 

For B = 2 

Un = 135.43 , U 33 = 86.59 and V n = 221.88 . (196) 

Comparing these numbers to those obtained in 7 using the exact numerical solution (Un = 127.8, 
U33 = 86.9 and Vu = 200.2), we see that the rational map ansatz has enabled us to obtain quite accurate 
moments of inertia. We recall that the old parameter set led to a model of the deuteron which was much 
too tightly bound. 

For B = 3 

u = 170.01 , v = 576.09 and to = -109.47. (197) 
These were evaluated in [5], using the exact numerical solution (u — 136, v — 435 and w = —91). 

For B = 4 

U n = 197.60 , C/33 = 236.49 and v = 911.45 . (198) 
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These numbers were calculated using the same procedure that we have used throughout, but with the 
old parameters. Walhout [28] performed a different style of analysis for the B — 4 Skyrmion, and 
unfortunately we are unable to directly compare our results for the individual components of the inertia 
tensors. 



C Appendix: Coefficients of Wavefunctions 

The FR constraints do not determine all coefficients in the wavefunctions. Usually finding these constants 
is trivial, but in some cases (as in the B = 5 first and second excited states) one has to be more careful. 
As an illustration let us consider the constants c± in (|94|) . The solutions of (191I92|I form a subspace of 
Hilbert space, which is transformed into itself when acted upon by the operator of the kinetic energy 
(fTTTj) . Therefore, the eigenvectors of the operator will define the wavefunctions we are looking for. In 
terms of the moments of inertia, c± is given by 



_ b 2 + 5&i - ai - 5a 2 ± -y/ (b 2 + 5bi - ai - 5a 2 ) 2 + 20(ai - a 2 )(&i - b 2 ) , . 

C± ~ 2V5(h-b 2 ) 



where 



1 / lOZ/u + 2V n + 20Wh 25C/ 33 + V 33 - WW; 



1 " 8{ U 11 V n -W? 1 + U 33 V 33 -W$ 3 )> ( ' UUj 
I ( 26U 11 + 2V 11 +4W 11 9U 33 + V 33 + 6W 33 \ 
G2 S\ U n V u - + U 33 V 33 - Wi 3 J ' 

1 / 10t/n + 2^11-4^11 25£/ 3 3 + U33- 10^33 \ f2()2l 

1 8\ U u Vn - U 33 V 33 - W| 3 ) ' ( ' 

. 1/ 26^1 + 2^1-20^! 9^33 + ^33 + 61^33 ^ 

2 8 ^ C/ U V U - W& U 33 V 33 - W 3 2 3 ) " 1 j 



The quantum energy of these states is given by 



E = ai+5a 2 + c ± \/5(bi-62) 
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